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Abstract: The entanglement entropy of an annulus is examined in a three-dimensional 
system with or without a gap. For a free massive scalar field theory, we numerically 
calculate the mutual information across an annulus. We also study the holographic 
mutual information in the CGLP background describing a gapped field theory. We 
discover four types of solutions as the minimal surfaces for the annulus and classify the 
phase diagrams by varying the inner and outer radii. In both cases, we find the mutual 
information satisfies the monotonicity dictated by the unitarity and decays exponen¬ 
tially fast as the gap scale is increased. We speculate this is a universal behavior in 
any gapped system. 
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1 Introduction 

Entanglement entropy is an invaluable tool to explore various aspects of quantum 
field theories (QFTs) in diverse dimensions, such as critical phenomena [1-8], confine¬ 
ment/deconfinement phase transition [9-13] and renormalization group flow [14-18]. It 
depends on a state of interest and the shape of an entangling surface E that divides a 
space into a region A and its complement A. In the simplest case, E is chosen to be 
a round sphere (or two endpoints of an interval in two dimensions) which allows us to 
use a conformal transformation [5, 19] and obtain analytic results for conformal field 
theories (CFTs). Multiple disjoint intervals for CFT 2 are examined by [20-24], and 
small deformations of an entangling surface in CFT^> 3 are perturbatively studied in 
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[25-28] recently. More general shapes, however, have not been fully understood so far 
because of the computational complexity, especially in non-conformal held theories . 1 

A special case is a system with a large mass gap where entanglement entropy can 
be expanded in powers of the inverse of the gap. The coefficients appearing in the 
expansion are unknown in general, but are assumed to be integrals of the functions 
of the extrinsic curvature and its derivatives of £ [31]. The entanglement entropy 
for £ cliffeomorphic to a circle is examined by [32] for free massive fields in three 
dimensions, showing that all the coefficients can be systematically determined by the 
logarithmic divergences of higher-dimensional theories that are the consequence of the 
conformal anomalies (see also [33]). Similar argument holds for £ diffeomorphic to 
multiple disjoint circles. 

In this paper, we consider the entanglement entropy of an annulus in QFT 3 as 
a guide to investigate the phase structure of the ground state . 2 For CFT, the finite 
part of the entropy is a function of the ratio of the inner and outer radii Ri < R 2 . 
Unfortunately, the conformal transformation used for a spherical entangling surface 
[19] does not help us identify the function for the annulus. It, however, was shown in 
[35] that the strong subadditivity [36] requires the function is concave with respect to 
\og(R 2 /R 1 ). In a gapped system, the finite part of the entropy depends not only on 
the ratio, but also on the gap scale, and thus there are no known constraints for the 
entropy from the strong subadditivity. 

To check if the constraints from the strong subadditivity holds for CFT 3 , or more 
generally to fix the dependence on the ratio and the gapped scale, we perform numerical 
calculations of the entanglement entropy for a free massive scalar held theory. We put 
the scalar held on radial lattice following [37] and compute the mutual information 
across the annulus (see Fig. 2). The mutual information is better than entanglement 
entropy itself in a sense that it is free from UV divergences and independent of the 
regularization scheme. For a massless scalar held, the strong subadditivity constrains 
the mutual information to be a convex function with respect to the ratio of the annulus. 
We conhrm that the convexity holds in our results and inspect the limits of R 2 /R\ —> 1 
and R 2 /R\ —>• 00 . In the former limit, we approximate the thin annulus by a thin strip 
(see Fig. 5) and evaluate the mutual information by dimensional reduction to an interval 
in (1 + 1) dimensions. In the latter case, we can conformally map the annulus to two 
disjoint circles whose entanglement entropy is studied both numerically and analytically 
in the large separation limit [38-40]. We fold that our fittings are consistent with the 
analytic results within our numerical precision. The implementation of the mass is 

1 See e.g. [29, 30] for studies on non-conformal theories including finite temperature cases. 

2 Refer to [34] as a related work on the fuzzy sphere. 
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straightforward numerically, and we observe that the mutual information exponentially 
decays as the mass increases while fixing the ratio R 2 /R 1 . 

Another model we are able to tackle is a strongly coupled QFT holographically dual 
to the gravity on the AdS space. The holographic calculation of entanglement entropy, 
known as the Ryu-Takayanagi formula [ 6 , 7, 41], associates the given region A in a 
QFT d to a codimension-two minimal surface 7,4 satisfying d'fA = FI in the AdS^+i space, 
and gives the entropy Sa by the area of the surface, Sa = Area( 7 u)/( 4 GV). It can be 
applied to disjoint regions and exhibits interesting transitions between different minimal 
surfaces with the same boundary condition [ 22 ], each one of them corresponding to a 
specific phase in the dual QFT. 

We study the holographic entanglement entropy of the annulus by extending the 
work [32] for a disk in a confining gauge theory with a gap described by the CGLP 
background [42]. The entropy given by the area of a minimal surface [ 6 , 7] shows a 
phase transition due to the change of the topology [17]. There are four types of minimal 
surfaces anchored on the annulus, namely, (1) hemi-torus, (2) two disk, (3) two cylinder, 
and (4) one disk and one cylinder (disk-cylinder) types. The last three solutions are 
superpositions of the disk- and cylinder-type solutions found in [32] . The first one was 
also constructed by [35] in the AdS space. Since the mutual information vanishes for 
disconnected surfaces only the hemi-torus solution has a non-zero value. Comparing 
their areas we classify the four phases in the (iR, i? 2 )-plane as shown in Fig. 11. The 
holographic model also exhibits the exponential decay of the mutual information with 
respect to the gap with the ratio R 2 /R 1 fixed. 

Based on the observations in the free massive scalar and the holographic models, 
we speculate that the mutual information through an annulus decays exponentially as 

^annulus ~ exp[-#m(R 2 ~ Rl)\ , (1-1) 

in any gapped system with a gap scale m. 

2 Entanglement entropy of annulus 

In this section, we review general properties of the entanglement entropy for an annulus 
A in CFT and a gapped system. We discuss its relation to the mutual information 
between the inner disk and the compliment of the outer disk. 
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2.1 Conformal field theory 

In a three-dimensional CFT, the form of the entropy for an annulus A is fixed by the 
conformal symmetry, 

S a (Ri,R 2 ) = a 2n{Rl + R2) - f{R 2 /Ri) , ( 2 . 1 ) 

where the first term obeys the area law with the UV cutoff length e and the second term 
/ is a function of the ratio R 2 /R\ of the radii. This function / should be monotonically 
decreasing and convex 


f\p) < 0 , f\p) > 0 , ( 2 . 2 ) 

with respect to the new variable p = \og(R 2 /R-i), due to the strong subadditivity 

Sb + Sc ~> Sbuc + Rbcc ■ (2-3) 

In what follows, we review the derivation of (2.2) given by [35]. 

Let the regions B and C be a disk of radius R 2 and an annulus of radii R\ and R 3 
with R\ < R 2 < Rj, as in Fig. 1 (a). The entanglement entropy for a disk of radius R 
takes a form of 

Sdisk(i?) = a 2 ^ - F , (2.4) 

with a constant F , 3 The strong subadditivity (2.3) together with (2.1) and (2.4) yields 
the monotonicity in ( 2 . 2 ): 


f(Rs/Ri) < f(R*/Ri) ■ (2.5) 

The convexity in (2.2) can be derived similarly by taking both B and C as an annulus 
of radii R 2 and /14. and an annulus of radii R\ and R 3 satisfying R\ < R 2 < R 3 < R 4 
as in Fig. 1 ( 6 ). In the 1 ? 4 — > R 3 limit, the strong subadditivity 

f(R 4 /R 2 ) + f(R 3 /Ri) < f(Ri/Ri) + }{Rz/R*) , (2.6) 

reduces to the monotonicity of 

3 The constant equals to the free energy on S 3 , F = — logZfS' 3 ) [19]. 
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(a) ( b ) 

Figure 1 . The subsystems B (in light gray) and C (in dark gray) to prove the monotonicity 
(a) and the convexity ( b ) of the function / in (2.1). The striped regions are the intersections 
BnC. 


2.2 A gapped system 

In theories with a mass gap of order m, the entanglement entropy of a region A has an 
expansion in powers of 1/m: 


Sa = ol -b /3rru 

e 


-7E + J] 


n =0 


L 2n+l 

m 2n+1 


(2.7) 


with numerical constants a, /3 and the topological entanglement entropy 7s [3, 4], Here 
the 7 s depends on only the topology of the entangling surface £ = dA and detects 
long-range order. The dimensionful coefficients cf/ +1 are postulated [31] as local inte¬ 
grals of functions of the extrinsic curvature and its derivatives on £. In other words, 
entanglement contributing to C 2 n+1 localizes on the entangling surface in the large-m 
limit due to the short correlation length of order 1/m. 

Applying (2.7) to the annulus A of our interest, the entanglement entropy should 
take the form of 

S a (Ri, R2, m) = a 27T<yRl + + 27 t/3 m (Ri + R2) ~ 7 s + V] ^n+i > ( 2 - 8 ) 

e L ' m + 
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where E is two concentric circles of radii R± and R 2 . The coefficients c% n+ i are polyno¬ 
mials of the radii of order —(2 n + 1). 

2.3 Mutual information 

The mutual information between two disjoint regions B and C is defined out of the 
entanglement entropies as 

I(B,C) = S b + S c -S B uC . (2.9) 

It is always finite because the area law divergences cancel by definition, and non¬ 
negative because of the subadditivity Sbuc < Sb + Sc- In addition, the strong subad¬ 
ditivity yields the monotonicity of the mutual information 

I(B, C) < I(B, C U D) , (2.10) 

for any region D. 

To extract the finite parts of the entanglement entropies (2.1) and (2.8) of the 
annulus A, we take B and C to be two regions outside A, namely, a disk of radius R\ 
and the complement of a disk of radius R 2 , respectively (see Fig. 2). We can interpret 
this mutual information as how much quantum information is shared by B and C across 
the annulus A. Since the entanglement entropy of a given region is equal to that of the 
complement, Sc and Sbuc equal the entropies of a disk of radius R 2 and an annulus 
of inner and outer radii R± and R 2l respectively. The mutual information / across the 
annulus A then reduces to 


I(R u R 2 ) = I(B,C) = Sdisk(fli) + S disk (R-2) ~ S a {R u Rv) . (2.11) 


In this setup, the inequality (2.10) translates into the monotonicity of / with respect 
to the radii Ri,R 2 \ 

R 2 )> 0 , R 2 ) < 0 , (2.12) 

which holds for any unitary QFT. The proof proceeds as follows: Let B,C be the 
regions in Fig. 2 and D be an annulus of radii R,\ + AR\ and R 2 , then the monotonicity 
I(B,C) < I(B U D,C) yields I(Ri,R 2 ) < I(R.\ + A Ri,R 2 ). Similarly let D be an 
annulus of radii R\ and R- 2 — A R, then the monotonicity I(B , C) < I(B, C U D) yields 
I{R\i R 2 ) A I(Ri, R-2 — AR 2 ). 

For CFT, the mutual information becomes 

IcFT = f(R 2 /Ri)-2F, (2.13) 
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Figure 2. The entangling regions for the mutual information. The region B is a disk of 
radius R \. The region C is the complement of a disk of radius R 2 . The complement of the 
union of the two regions B U C is the annulus A in the red colored region. 

with the constant F in the disk entropy (2.4), and the inequalities (2.12) are equivalent 
to the monotonicity of / that was already derived in ( 2 . 2 ). 

On the other hand, applying (2.8) for a gapped system to (2.11) leads to 


(2.14) 


-^gapped TS 2y r ]i sk • 


Here the m-dependent terms cancel out due to the assumption that the dimensionful 
coefficients cf n+1 in (2.8) are integrals on the entangling region £. Note that the 
expression (2.14) would fail for small masses such as mR\ < 1 or m(i ?2 — Ri) 1$ 1 if 
there could exist an exponential term like O (exp [—#m]) to (2.7) which can not be seen 
in the large mass expansion. We will discuss such a correction in Section 5. 

In the following sections, we will use these mutual informations (2.13) and (2.14) 
to determine the function / in CFT and to check whether the large mass expansion 
formula (2.7) holds for the annulus. 

3 Free massive scalar field 

Let us apply the general discussion on the annulus entropy in Section 2 to a free massive 
scalar held whose action is defined by 



(3.1) 
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In this case, the coefficients f3 and 7 in the entropy (2.8) are known to be f3 = —1/12 
and 7 = 0 . 1 The coefficients C 2 n+1 are calculated [32, 33] up to n = 1, being local 
integrals of functions of the extrinsic curvature k and Fs derivatives on the E. For 
example, 


Cl 


n 0 + 3ni /2 
480 



(3.2) 


for no free scalar fields and rii /2 free Dirac fermions. I 11 the present case, the en¬ 
tangling surface is two disjoint disks of radii R±, R 2 whose extrinsic curvatures are 
k = 1 /Ri, l/i? 2 - Thus cf = — 2 io(l/-Ri + I/-R 2 ) f° r a single free scalar field. The con¬ 
stant term F of the disk entropy (2.4) is analytically calculated as the free energy on 
a three-sphere, F sc a i ar = (ln2)/8 — 3C(3)/167 t 2 ~ 0.0638 [43]. The mutual informations 
(2.13) and (2.14) are 


‘massless 


L massive 


/ (R2/Ri) — 2 F sca \ ar , 

0 . 


(3.3) 

(3.4) 


3.1 Numerical results 

We perform the numerical calculation by putting a free scalar field on the radial lattice 
following [37, 44], whose details can be found in Appendix A. The main results are 
presented in Fig. 3 and 4. 

Fig. 3 shows the mutual information / (3.3) for the free massless scalar field. The 
function / satisfies the desired monotonicity and convexity ( 2 . 2 ) with respect to p = 
log{R 2 /Ri) as is clear in Fig. 3 (a). The mutual information / asymptotically vanishes 
for large R 2 /Ri, which means that the function / has a finite constant term 2F sca i ar . 
This suggests that the finite constant term is topological and additive for each connected 
component of the entangling surfaces, namely, proportional to the 0-th Betti number 
60 [£] ofS. The numerical function / = J(i? 2 /Ri) is well approximated by h/(R2/R\ — l) 
with h ~ 1/4 for small width, but h monotonically increases to h ~ 1/3 for large R2/R1 
(see Fig. 3 ( 6 )). We therefore propose that / is given by 

/(ft/ft) = /WTA + 2, (3.5) 

where h(R 2 /Ri) is a mild monotonically increasing function of R 2 /Ri such that /r ~ 1/4 
for R 2 /R\ ~ 1 and h ~ 1/3 for R 2 /R\ 1. These asymptotic values are consistent 

with previous works [40, 45] as will be explained in the next subsections. 

4 The topological entanglement entropy vanishes because there is an empty theory in the IR of the 
massive scalar theory. 
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(a) (&) 

Figure 3. The mutual informations I across the annulus of radii R± and R2 for the free mass¬ 
less scalar field, (a) The mutual information I (the orange line) has the desired monotonicity 
and convexity, and is well fitted by h/(R, 2 /R\ — 1) (the black dotted line). ( b ) However, this 
coefficient h is not a constant and increases with R2/R1 from h ~ 1/4 (the black dotted line) 
to h ~ 1/3 (the blue dotted line). 

The mutual information (3.3) for the free massive scalar field is displayed in Fig. 4. 
It is monotonically decreasing with the mass (i.e., decreasing with mi? 2 or mR\ while 
R2/R1 being fixed), and almost vanishes for large mass (Fig. 4 (a)) as is consistent with 
(3.4). In fact, Fig. 4 (6) demonstrates that the mutual information decays exponentially 
with a “dimensionless width” m(i ?2 — Ri), 

/massive oc m(R 2 + aRi) exp [-bm(R 2 - Ri)} , (3.6) 

with constants a and b. This exponential behaviour satisfies the expected monotonicity 
(2.12). We will find similar decay even in the holographic model in Section 4 and discuss 
their possible universality in a gapped phase in Section 5. 

3.2 Small and large width limits in CFT 

The annulus with small width {R 2 /R\ ~ 1) can be approximated by a thin strip of 
width i ?2 — R\ extending along a circle of radius 2 tiR\ as in Fig. 5. 5 The mutual 
information for the thin strip of width 5 is shown to obey [45, 46] 

1 ~ K J ’ ( 3 - 7 ) 

5 We thank T. Takayanagi for drawing our attention to this point. 


9 
















Figure 4. The mutual informations I across the annulus of radii R\ and R -2 for 
scalar fields with different masses m. (a) I monotonically decreases with the mass m 
(orange—ffilue—)>green—»red). ( b ) In fact, it exponentially decreases with the dimensionless 
width m(i?2 — Ri). For ?n(i?2 — Ri) > 1, it shows I oc m(R .2 + aR\) exp[—6m(i?2 — Ri)] with 
a~2~5(a = 3in the figure) and b ~ 2.5. 


where A is the area of the plane bounding the strip. This behavior was derived by 
dimensionally reducing the thin strip to an interval in (1 + 1) dimensions for free fields 
and summing the mutual informations over the Kaluza-Klcin modes. 

The coefficient k is calculated for a free massless scalar held [45] to be k = 0.0397. 
Applying (3.7) to our case, we find 


/ ~ 0.0397 


2vrRi 


0.249 


R 2 - A i Ri/Ri-l’ !>N; 

which fits our numerical result in the small width limit (h ~ 1/4 in (3.5)) very well. 
One may wonder if the small width limit of the mutual information (3.7) is universal 
and the coefficient k counts the number of degrees of freedom in any QFT. We will 
come back to this point in Section 5 where we calculate k in a holographic model. 

Next, consider the opposite limit where the width is large. Let z t (i = 1, 2) be 
the two-dimensional Cartesian coordinates related by an inversion transformation 


(+ %o)i — Rt 


2 (w - Wo)i 


\W — Wo 


(3,9) 


where Wq is the inversion point. The inverse map is obtained by exchanging the role 
of w and z in the transformation with the inversion point at z = z$. Rt is a constant 
which we can tune arbitrarily. 

Consider an annulus in the w-coordinates whose center is at the origin with radii 
Ri < R 2 - Let the points at W 2 = 0 on the outer circle be Pi, p 2 and on the inner circle be 
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(a) 


Ri 


R2 — Ri 

(b) 


Figure 5. A thin annulus (a) can be approximated by a thin strip (b) with compactified 
direction. 


qi,q 2 ■ We choose the inversion points w$ and Zo on the real axes at (w\,w?) = (i?o,0) 
and (z\, Z 2 ) = ( R' 0 ,0), respectively. We assume wo is inside the annulus, R± < Ro < Ro¬ 
under the transformation (3.9), the annulus is mapped to two disjoint circles 6 7 (see 
Fig. 6) and the points p \, p 2 and q \, g 2 are at the intersections of the real axis and 
circles of radii R\ and R’ 2 given by 


R[ = R 2 7 


Ri 


Rl 


Ri 


2 > 


RL = Rl 


R2 


Ri 


Rl 


The distance between the centers of the two circles is 

Ri 


r = R^Ro 


R\ 


(Rl - Rl)(Rl - R\) ' 


(3.10) 


(3.11) 


The conformal symmetry implies that the cross ratio' x is invariant under the conformal 
transformation, 


\P1-P2Wq1 ~ gg] 

\pi ~ q2\\p2 ~ qi\ ' 

6 We thank K. Ohmori and Y. Tachikawa for the discussions on this map. 

7 There are two cross ratios for four points. The other one is 

= \pi ~P 2 |ki ^ 92! 

\pi - Qi\\P 2 - 92 1 


(3.12) 
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which in our case is 


4R X R 2 _ 4 R[R' 2 

(Ri + R 2) 2 ~ r' 2 - {R[ - R’ 2 f ' 


(3.13) 



Figure 6. The inversion map of an annulus to two disjoint circles. The region inside the 
annulus in red color is mapped to the outside of the two circles in red color. 


In this way, we can calculate the entanglement entropy of two disjoint circles from 
that of the corresponding annulus. The former was studied in [38-40] in the widely 
separated limit for a free massless scalar held. The mutual information between the 
two circles is [40] 

1 = \^T + 0 (WK/’-' 2 ) 2 ) . (3.14) 

The inversion maps (3.10) and (3.11) convert it to the mutual information of the an¬ 
nulus, 


/ 


1 1 

3 rJWR'" ’ 


(3.15) 


in the large width limit (R 2 /R\ 1). What we observed in the previous subsection is 

nothing but this asymptotic form consistent with numerical result shown in Fig. 3. 


4 Holographic entanglement entropy 


In this section, we examine the entanglement entropy of an annulus in CFT 3 and a 
gapped system holographically described by the Einstein-Hilbert gravity in the (asymp¬ 
totically) AdS 4 space. The holographic formula [6, 7] 


Sa = 


min 


Area [7 a] 
4 Gn 


7 


(4.1) 


- 12 - 





















associates the entropy of a given region A to the area of a codimension-two minimal 
surface 7 a homologous to the region A, i.e., djA = £• Fig. 9 illustrates the cases for A 
being a disk. 

If there are multiple extremal surfaces, we always pick one of them with least area 
according to the formula (4.1), which yields a transition between minimal surfaces as 
we vary a parameter such as a gap scale. In this sense, each extremal surface can be 
regarded as a phase in QFT as we will see in the following. 


4.1 The AdS 4 background 

We start with CFT 3 dual to the AdS 4 background 


2 9 dz 2 — dt 2 + dr 2 + r 2 dd 2 

ds z = L --- 


(4.2) 


with the AdS radius L. The original CFT 3 is interpreted to live on the boundary z = 0 
(or at z = e <C 1 if UV regularization is needed). 

The extremal surface respecting the rotational symmetry of the annulus is a solution 
to the equation of motion for the action 


I[r(z)\ 


7 tL 2 


2 G 


N 


A r{z)yj 1 + r'(z ) 2 
* ’ 


(4.3) 


with the boundary conditions r(0) = Ri [i = 1,2) on its ends. There are two possible 
extremal surfaces depending on their topologies: 


• Two disk phase (Fig. 7 (2)): 7^4 is the superposition of disconnected two disks, 
each of them being given by 

r (z) = sjR 2 - z 2 , (i = 1,2) , (4.4) 

respectively. This solution always exists independent of the size of the annulus. 

• Hemi-torus phase (Fig. 7 (1)): 7,4 is a connected extremal surface. The analytic 
solution is obtained in the following way [47-49]. It consists of two branches in 
the (r, z)-plane as 

^ fhhexp [~f-(z/r)\ , 

[R 2 exp[-f + (z/r)} , 
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where the functions f±{x) are defined using the incomplete elliptic integrals 8 by 


f±(x) = ^log(l +x 2 ) ±rjx m [F (uj(x)\t) 2 ) - II (l - r) 2 ,u(x)\r) 2 )] , (4.7) 


with the range 0 < x < x m 



= arcsin 


_ x/Xm _ 

y/l—r) 2 (l—x/x m ) 


The 


parameter 77 in the range 77 e [1 /a/ 2, 1] is related to the ratio R 2 /R 1 of the inner 
and outer radii of the annulus as 


\og(R 2 /Ri) = 277 


2?7 2 - 1 
1 — 77 2 


[KO 2 ) - n(i - >) 2 |i) 2 )] 


This solution is available only for (1 <)R 2 /R\ < 2.724. 


(4.8) 


The two disk phase is realized for the large width R 2 /R\ > 2.724 where it is the 
unique solution, while it compete with the hemi-torus phase when R 0 /R 1 < 2.724. In 
order to fix the location of the phase transition, we calculate the mutual information 
/ across the annulus defined by (2.11). It is clear in the holographic setup that I > 0 
signifies the hemi-torus phase because / = 0 in the two disk phase. 9 We benefit from 
the relevant result of [49] to get the mutual information in the hemi-torus phase 10 


7 tL 2 


*■ hemi-torus 


G 


N 


E(r] 2 ) - (1 - 77 2 )K(t7 2 ) 
y/ 2 V 2 -1-1 


(4.10) 


whose plot is displayed in orange color in Fig. 8. It is a two-valued function whose 
lower branch is always negative and the upper branch intersects with / = 0 at R 2 /R\ = 
(R 2 /Ri) critical ~ 2.4. Since the holographic formula (4.1) selects the non-negative /, the 

s The definitions of the incomplete elliptic integrals used here are 

fX 1 

F(a;|m) = / d8 — — , 

7o \J\ — m sin 2 0 

nx x ( 4 - 6 ) 

II(n, x\m) = / dd - ^- , 

Jo (1 — nsin 2 8) \J 1 — msin 2 6 


and K(m) = F(7 t/2|to) and II(u|77i) = II(n,7r/2|m). 

9 The mutual information can vanish only in the large -N limit and there are 0(1/N) corrections 
[50, 51] for finite N. More generally, the mutual information is bounded from below. 

10 The elliptic integral of the second kind is defined by 


E(m) 



(4.9) 
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(1) Hemi-torus phase 


(2) Two disk phase 


Figure 7. Two phases for the minimal surface in the AdS 4 background: connected hemi-torus 
phase (1) and disconnected two disk phase (2). Here the time t direction is suppressed. 


physical mutual information is given by /hemi-torus for R2/R1 < {R2/Ri) critical and 1 = 0 
for (R2/Ri)critical < R2/Ri- It has a kink at R2/R1 = (/^/-Ri) critical caused by the 
phase transition of the extremal surface 7 , 4 . Comparing with the general form (2.13) of 
the mutual information in CFT, Fig. 8 demonstrates the monotonicity and convexity 
( 2 . 2 ) of the function / with respect to p = log(R 2 /Ri). In other words, the holographic 
entanglement entropy of an annulus satisfies the strong subadditivity as guaranteed by 
the holographic proof based on the minimality of the surfaces [52] . 

4.2 The CGLP background 

We move onto a gapped theory described by an asymptotically AdS geometry whose 
IR region (away from the boundary) is capped off. As a concrete example, we use the 
CGLP background [42] in M-theory dual to a (2 + l)-dimensional QFT with a gap 
scale. 

The CGLP background is a (3 + l)-dimensional geometry times a seven-dimensional 
internal manifold, which asymptotes to the AdS 4 space times the Stiefel manifold V 5)2 - 
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log R 2 /R] 


Figure 8. The holographic mutual information I across the annulus of radii R± and R 2 for 
CFT. This mutual information I = /(R 2 /A 1 ) has a phase transition at (-R 2 /-R 1 Critical — 2.4, 
and vanishes for R 2 /R 1 > (R 2 / Ri) critical because the disconnected two disk phase is realized. 


In the Einstein frame, the metric is given by 


ds 2 = a(u) [du 2 + /3(u) (— dt 2 + dr 2 + r 2 d6 2 )~\ + gijdy l dy j , (4.11) 


where u is the holographic coordinate of the AdS 4 ranging from the IR capped-off point 
0 to the UV fixed point 00 . y l (i — 1 , - • • ,7) are the coordinates of the internal manifold 
with a volume 

V{u) = j Y[dy l t/detg , (4.12) 

J i=1 


vanishing at u — 0. The functions in the metric are given by 


a(u) 

V(u) 

H(n) 



,, 4 

4 ’ H(u)c 2 (u) 

—-——— H‘^{u)(2 + cosh u) 3 / 8 sinh 3,/2 j sinh 3 ^ 2 u , 

~^ 6 q3/2q11/4 f & 

e ^ 2 ./(2+coshu) 1 / 4 ( t 4 - l) 5 / 2 ’ 

37/ 4 £ 3/2 (.QgpS (u/2) 

2(2 + cosh u) 3 / 4 


(4.13) 


with two dimensionful parameters L and e. The parameter L is the AdS radius near 
the boundary, determined by the number of M2-branes N and the Planck length i v as 
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L = 3 _2 / 3 27T 1 / 3 £ p A rl//6 . The parameter e, defining the size of deformation [53], has mass 
dimension —4/3, letting H be dimensionless. V appears to depend on e, but does not 
indeed. By rescaling the boundary coordinates (t, r ) appropriately, one can remove e 
completely from the metric if one wishes. 

Let us take a look at the UV behavior of the metric (4.11) for a moment. When u 
is close to the UV cutoff m 4 A > 1, the function H(u) becomes 

H{u) ->• 2 15/4 3 3/4 L 6 e _9 “ /4 , (4.14) 


and the other functions approach 

a(u) -> — L 2 , j8(u) -> 2 3/2 3“ 5/2 L “ 6 e 3 “ /2 , 

16 

V(«) -4 3VL !1 ' ! , c 2 ( m ) -> 2- 13 ' 4 3 7 / 4 e s “/ 4 . 


(4.15) 


The transformation z = 2 5 / 4 3 1,/4 L 3 e 3 “/ 4 takes the metric to the Poincare coordinates 
of the AdS 4 space near the boundary 


2 r2 dz 2 — dt 2 + dr 2 + r 2 d # 2 

ds 2 L 2 ---1-. 

z 1 


(4.16) 


Since the extremal surface for a small annulus localizes near the boundary, the 
entanglement entropy remains to have the previous two phases shown in Fig. 7 in the 
CGLP background. In addition, there are new phases for a large annulus whose minimal 
surfaces can reach and terminate on the IR cap-off as we describe below. These are 
superpositions of disk- and cylinder-type solutions for a disk region [17, 32] depicted 
in Fig. 9. They have different topologies as the names suggest, and the cylinder-type 
solution only exists and dominates for a large radius. This resembles the situation for a 
strip region in a gapped system, which is interpreted as a confinement/deconfinement 
phase transition [9, 10]. In the present case, the minimal surface switches from the disk- 
type to the cylinder-type at the critical radius R = ^critical — 0.72/m, where m = e~ 3 ^ 4 
is the gap scale determined by the CGLP metric. Taking into account these facts, 
we end up with three superposed phases; two disk phase, one disk and one cylinder 
(disk-cylinder) phase, and two cylinder phase. The first one has already appeared for 
CFT in the previous subsection (see Fig. 7). The second and third ones are drawn in 
Fig. 10. In total, there are the four phases for the annulus in the CGLP background: 


(1) the hcmi-torus phase (Fig. 7 (1)) for R 2 — R± < 1/m. 

(2) the two disk phase (Fig. 7 (2)) for Ri, R 2 < ^critical, 

(3) the disk-cylinder phase (Fig. 10 (3)) for Ri < R cv i tica i < R 2 , 
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(a) disk phase 


(b) cylinder phase 


Figure 9. Two phases of the extremal surface in calculating holographic entanglement en¬ 
tropy of disks in the CGLP background: disk phase (a) and cylinder phase ( b ). In the 
Poincare coordinate z = 2 5 / 4 3 1 / 4 L 3 e~ 3 “ //4 , the UV boundary u = A corresponds to z = e = 
2 5 / 4 3 1 / 4 T 3 e -3A / 4 and the IR capped-off point u = 0 corresponds to z = zq = 2 5 / 4 3 1 / 4 h 3 . In 
the cylinder phase, the extremal surface terminates on the IR capped-off point z = zq. 


(4) the two cylinder phase (Fig. 10 (4)) for i? critica i < Ru R 2 , 

There are apparently overlaps between the first phase and the others, where the 
one with the least entropy is realized. To classify the phase structure, we calculate the 
holographic entanglement entropy in a similar way to the previous pure AdS case. The 
rotational symmetry lets us assume the radial coordinate r of the extremal surface ja 
as a (two-branched) function r± = r±(u ) of the holographic coordinate u. The area 
functional becomes 

J [ r ( M )] = 2^^ / dur ±( u ) + P(v)(r±(u)) 2 , (4.17) 

with < 7 (it) = V(u)a(u)/3 1 / 2 (u). The extremal surface r = r(u ) should satisfy the equa- 
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(3) Disk-cylinder phase 


(4) Two cylinder phase 


Figure 10. Two new disconnected phases for the minimal surface in the CGLP background: 
disk-and-cylinder phase (3) and two cylinders phase (4). In the Poincare coordinate 2 = 
2 5 / 4 3 1 / 4 L 3 e _3n / 4 , the UV boundary u = A corresponds to z = e = 2 5 / 4 3 1 / 4 T 3 e~ 3A / 4 and the 
IR capped-off point u = 0 corresponds to z = zq = 2 5 / 4 3 1//4 L 3 . 


tion of motion 


2g(u)y/l + f3(u)(r'(u)) 2 


r(u)g(u)/3(u)r'(u) 

V 1 + P ( u )( r ’( u )) 2 


(4.18) 


with the boundary conditions r + (oo) = R 2 and r_(oo) = R 1 . In contrast to the CFT 
case, the analytic solution remains to be known. Instead, we employ the numerical 
calculation by the “shooting method”. 


• In the hemi-torus phase, we solve the equation of motion (4.18) from the tip 
(r,u) = where the two branches meet and have an expansion 


f±{u) = r* ± 2 \ 


g(u* 


g(u*)/3'(u*) + 2 g'(u*)/3(u*) 


\/u — u* + 0({u — m*) 3 / 2 ) . (4.19) 


The radii of the annulus (R 2 , Ri) = ( r + (u = oo),r_(w = oo)) are functions of 
(r*,n*), respectively. 
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• In the three disconnected phases, the extremal surfaces 7 a for the annulus are 
obtained just by summing the two extremal surfaces 7 disk(_Ri) and 7 disk(j? 2 ) f° r two 
disks of radii R± and R 2 . The extremal surface 7 disk(/?) for a disk in the CGLP 
metric was obtained [32] as follows. The disk-type solution can be constructed 
by solving the equation of motion (4.18) from the tip of the disk (r,u) = (0,w*), 
where the extremal surface shrinks as 


r(u) = 


2 g(u* 


2 /3(u*)g'(u*) + g(w*)/3(w* 


-a Ju — u* + 0((u — u*f 2 ) . (4.20) 


On the other hand, the cylinder-type solution extends to the IR capped-off point 
u = 0 and we solve the equation of motion (4.18) from (r,u) = (r*,0) where the 
extremal surface terminates and behaves as 

r(u) = r * + » 7 ( o )“ 2 + 0( “ 3) ' (4 ' 21) 

The disk radius R = r( 00 ) is given as a function of u* or r*, respectively. 


After solving the equation of motion numerically, we compare the holographic 
entanglement entropies (4.17) between the four phases. The resulting phase diagram is 
presented in Fig. 11. It shows that the hcmi-torus phase is realized when the width of 
the annulus is small against the gap scale. Note that there is no phase for R 2 /R\ < 1 
since R 2 is the outer radius of the annulus. 

The mutual information (2.11) across the annulus vanishes in all the disconnected 
phases, and I > 0 only in the hemi-torus phase. Fig. 12 shows / as a function of 
log(i? 2 / Ri) with mR 2 fixed. It is positive and decreases as R 2 /R\ becomes large, 
but vanishes at some point due to the phase transition from the hemi-torus phase to 
a disconnected phase. It is also monotonically decreasing with the mass for a fixed 
R 2 /R\. We will discuss the mass dependence of the mutual information in the next 
section. 


5 Universal behaviors 

In the last two sections, we have dealt with the annulus entropies Sa{Ri, R 2 ) or the 
mutual informations / across the annulus for the free massive scalar theory and the 
holographic model. In this section, we will compare these two cases, and attempt to 
identify universal behaviors of entanglement entropy. 

First we consider the small width limit of the mutual information in CFT. From 
the held theory result, we anticipate (3.7) holds even in the holographic model. Since 
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mR i 


Figure 11. The phase diagram of the entanglement entropy for an annulus of radii R± and 
i? 2 - The hemi-torus phase is favored when the width of the annulus is small compared to the 
gap scale. 
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Figure 12. The holographic mutual information / across the annulus of radii R\ and R 2 in 
a gapped theory holographically described by CGLP metric. I = I(rnR\. rnR/ 2 ) vanishes for 
large R 2 /R 1 because the phase becomes disconnected, and monotonically decreases with the 
mass m increased. 

the hemi-torus phase is always favored in the small width limit, we can make use of 
the relations (4.8) and (4.10). A short calculation yields the small width behavior (3.7) 
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with the coefficient K^oi given by 11 


_ L 2 r[ 3 / 4] 4 
Kho1 = 2 xG n 


(5.1) 


It is plausible that k in (3.7) counts the effective degrees of freedom in a given QFT 
because it is proportional to the number of fields in free held theories which characterize 
the UV fixed point detected by the small width limit of the mutual information. Indeed, 
the Khoi hr the holographic model decreases under any RG how thanks to the holographic 
c-theorem [15, 16, 55, 56] that provides the constraint Luv > Air for the AdS radii 
in the UV and IR hxed points. Similar story may hold for the mutual information 
through two concentric (d — 2)-sphere separated by a short distance 6 which behaves as 
I ~ /t Area (S d ~ 2 )/S d ~ 2 in d > 4 dimensions [45, 46]. We do not explore this possibility 
in this paper, but hope to investigate it in the future. 

In a gapped system, we observed the exponential decay of the mutual information 
(3.6) for a free massive scalar held. It provides a strong evidence for the validity of the 
ansatz (2.7) of the entanglement entropy expanded with respect to the inverse of the gap 
scale, whose coefficients are the integrals of local invariants localized on the entangling 
surface. It also implies the existence of an exponentially suppressed correction to the 
ansatz that will never be seen in the large gap expansion. It is of interest to see to 
what extent the ansatz (2.7) captures the feature of entanglement entropy in a gapped 
system. Actually, our holographic calculation in the CGLP background exhibits the 
exponential decay of the mutual information as in Fig. 13. These observations suggest 



-2 


mR2 

+ 1.0 

X 1.5 



m(R 2 -Ri) 


Figure 13. The exponential decay of the holographic mutual information I with the dimen¬ 
sionless width m(R2—R±). For 0.1 < m(R-2—Ri) < 0.3, it shows I oc mR 2 exp[— b' m(R2~ 7?i)] 
with b' ~ 10. For m(i ?2 ~ Ri) > 0-3, this exponential behavior ends because of the phase 
transition to the disconnected phases with 1 = 0. 

11 See [54] as a recent related work. 
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that the entanglement entropy in a system with a gap m has a power series expansion 
of 1 /m with an exponential correction 

b 00 c s 

S A = - 7s + ^2 m^n+i + 0 (exp[-m8]) , (5.2) 

n =0 

where 5 is proportional to the shortest distance between disjoint entangling surfaces. 
This is equivalent to the speculation (1.1) for the mutual information of the annulus 
where 8 oc R 2 — Ri- We conjecture that (5.2) is a universal property in any gapped 
system. This resembles the universal thermal corrections in entanglement entropy [57- 
61] and it would be intriguing to find a relationship between them. 
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A Details of numerical calculations 


In this appendix, we summarize the numerical algorithm for calculating the entangle¬ 
ment entropy of the annulus for a free massive scalar held whose action is given by 
(3,1). 

A.l Radial lattice discretization 

We use the polar coordinates to put the theory on the radial lattice 

ds 2 = — dt 2 + dr 2 + r 2 dd 2 . (A.l) 


The radial coordinate r is discretized to N points with lattice spacing a. After the 
Fourier decomposition along the angular coordinate 9, the lattice Hamiltonian becomes 


H 


1 

2 


X 


' N 

X 

i =1 




N 

X 

l,j =1 




n,j 


(A.2) 


where and 7T n ., are the discretized scalar held with angular momentum n on the 
i-th site and its conjugate, respectively. The matrices K// depend on the angular 
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momentum and the mass m 

Q n 2 

K ’ 1 = - + n 2 + (ma ) 2 , K? = 2 + - + (ma) 2 , 

2 

These are related to the two-point functions of the scalar fields (X n ) t j = (4> n ^4> n ,j) and 
the momenta (P n )ij = (^n,i^n,j) as X n = \Kn 1/2 and P n = \Kl/ 2 . 

The outer and inner radii of the annulus are chosen to be half-integers in units of 
the lattice spacing, Ri/a = ri + 1/2 and i? 2 /a = r 2 + 1/2 with integers r 1; r 2 . This choice 
corresponds to the free boundary condition in the continuum limit. In our calculation, 
we vary r 2 from 100 to 120 and r± from 5 to r 2 — 5. The entanglement entropy of the 
annulus S(R\, f? 2 ) is obtained by using (r 2 — ri) x (r 2 — ri) submatrices (X^ 1,r2 )ij and 
(Pn 1,T2 )ij of the correlation functions X n , P n with the ranges r 1 + 1 < i, j < r 2 as 

OO 

S(R ll R 2 ) = S 0 + 2j2S n , (A.4) 

n =1 

where S n is the contribution from the n-tli angular mode 

S n = tr [{C n + 1/2) \og{C n + 1/2) - {C n - 1/2) log(C7 n - 1/2)] , (A.5) 

with C n = a /Xn 1 ,r 2 Pn 1,r2 ■ In the following, we describe how to perform this infinite 
summation over n under controlled numerical errors. 

A.2 Finite lattice size effect 

To avoid the finite lattice size effect, we repeat the calculation of S n (A.5) by changing 
the lattice size N and fit the results S n (N ) with the asymptotic expansion for large N 

^max 

S„(iV) = S n ( oo) +Y.Jfk- ( A '6) 

k =1 1 

We then read off the constant part S^oo) as the value of S n in the large-fV limit. 
Starting from N = 200, we increase the lattice size by AN = 20 until the resultant 
S n (oo ) stops changing up to error 5 = 10 -6 . We choose the fitting parameter /c max so 
that the maximum lattice size N is as small as possible. Typically we find A: niax = 3 ~ 
10 . 

The finite lattice size effect dominates only for small angular momenta n with small 
masses ma. In our calculation, the maximum lattice size reaches N ~ 1000 for n < 10 
in the massless case, but N = 200 is sufficiently large for n > 20 or ma > 0.1. The 
total numerical error in (A.4) can be estimated to be 0(206) < O(10 -4 ). 


k: 


1-1 


= K i+1,i = 


i + 1/2 

VW+ 1 ) 

(A.3) 
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A.3 Large angular momentum 


In the large angular momentum limit n —> oo, the correlation matrices X n and P n 
approach almost diagonal matrices [32]. The products of the submatrices X^ 1,r2 P^ 1,r2 
almost equal to 1/4 times unit matrix up to order 1/n 8 . The nontrivial entries are at 
the upper-left corners 


^ ,r 2 pr\,r 2 \ri+l,ri+l 

( ^ri,r 2 pri,r 2 ^r\+l,r\+2 


1 rl(ri + l) 2 rf(r 1 + l) 2 (2r! + l) 2 (m 2 + 2) 


4 16n 4 

r 3 (ri + l) 3,/2 (ri + 2) 3 / 2 
64n 6 

r 2 (ri + 2) 3//2 (ri — l) 1</2 (3(ri + l) 2 — 1) 
64n 6 


32 n 6 
+ 0(l/n 8 ) , 


+ 0(l/n 8 ) , 


+ 0(l/n 8 ) . 


(A.7) 


and at the lower-right corners 

( X r„r 2p r„r 2 y 2 ,r 2 _ 1 4(r 2 + l) 2 _ r 2 (r 2 + l) 2 (2r 2 + l) 2 (m 2 + 2) 8 

1 n n ’ “4 + 16n 4 32n 6 + [/ 


(X^ipri^^-l = r 2 (^2 ± ^ + 0(1/n 8 

r 2 /2 (r 2 - l) 3/2 (^2 + l) 2 (3r| - 1) 


W’ r2 C r2 ) 


,r 2 \r 2 -\,r 2 _ 


64 n 6 


+ 0(l/n 8 ) . 


(A.8) 


Here we restrict the ranges of ri,r 2 to 3 < ri and r\ + 3 < r 2 to avoid the overlap 
between the upper-left and lower-right corners, which is satisfied in our set up with 
5 < rq < r 2 — 5. 

The r 2 — rq — 2 eigenvalues of the matrix y/Xp 42 Pn 1 ,r2 are 1/2 + 0(l/n 8 ) and the 
other two are given by 


1 + >) _ (q) ( 2r ° + 1 ) 2 ( ? » 2 + 2 ) 

2 + n n 2n 2 


(a) = r a(^a + 

n - 16n 4 


a = 1,2. (A.9) 


Therefore, most of the eigenvalues do not contribute to the n-th entanglement entropy 
(A.5) up to order 1/n 8 and we obtain 


s»= E 


0 = 1,2 


•>(“) 


(y i ( a)\ (2r a + l) 2 (m 2 + 2) / \ 

(1 - log 4 a) ) + -- XX - -cx> log c. 


2 n 2 


(a) 

n 


+ 0(l/n 8 ). (A.10) 


This asymptotic formula is much faster than the direct calculation of (A. 5). 

We perform the matrix trace calculation (A. 5) for n less than some large angular 
momentum n*, and use this asymptotic formula (A. 10) for n > n* as long as S n (— 
0(logn/n 4 )) is larger than the machine precision. The other higher modes are ignored. 
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Our n* is determined as follows. Let the error of 0(l/n 8 ) in (A.10) be /z/n 8 
with n = fj,(m,ri,r 2 ). Then the total numerical error in (A.4) is estimated to be 
£“(/V™ 8 ) n/(7nl). We take n* to be the angular momentum where the asymp¬ 

totic formula (A. 10) agrees with the matrix trace calculation (A.5) up 7 8/n. Then 
/Lt/rz® < 75/n* holds and the total numerical error in (A.4) is bounded by X^(/V n8 ) ~ 
fi/(7nl) < 5. In this way, we can handle the numerical error within 0(6). 

References 

[1] C. Holzhey, F. Larsen, and F. Wilczek, Geometric and Renormalized Entropy in 
Conformal Field Theory , Nucl.Phys. B424 (1994) 443-467, [hep-th/9403108], 

[2] G. Vidal, J. Latorre, E. Rico, and A. Kitaev, Entanglement in Quantum Critical 
Phenomena, Phys.Rev.Lett. 90 (2003) 227902, [quant-ph/0211074]. 

[3] A. Kitaev and J. Preskill, Topological Entanglement Entropy, Phys.Rev.Lett. 96 (2006) 
110404, [hep-th/0510092], 

[4] M. Levin and X.-G. Wen, Detecting Topological Order in a Ground State Wave 
Function, Phys.Rev.Lett. 96 (2006) 110405. 

[5] P. Calabrese and J. L. Cardy, Entanglement Entropy and Quantum Field Theory , J. 
Stat. Mech. 0406 (2004) P06002, [hep-th/0405152], 

[6] S. Ryu and T. Takayanagi, Holographic Derivation of Entanglement Entropy from 
AdS/CFT, Phys.Rev.Lett. 96 (2006) 181602, [hep-th/0603001], 

[7] S. Ryu and T. Takayanagi, Aspects of Holographic Entanglement Entropy, JHEP 0608 
(2006) 045, [hep-th/0605073] . 

[8] T. Grover, A. M. Turner, and A. Vishwanath, Entanglement Entropy of Gapped Phases 
and Topological Order in Three Dimensions, Phys.Rev. B84 (2011) 195120, 

[arXiv: 1108.4038], 

[9] T. Nishioka and T. Takayanagi, AdS Bubbles, Entropy and Closed String Tachyons , 
JHEP 0701 (2007) 090, [hep-th/0611035], 

[10] I. R. Klebanov, D. Kutasov, and A. Murugan, Entanglement as a Probe of 
Confinement, Nucl. Phys. B796 (2008) 274-293, [arXiv: 0709.2140]. 

[11] A. Pakman and A. Parnachev, Topological Entanglement Entropy and Holography, 
JHEP 07 (2008) 097, [arXiv: 0805.1891], 

[12] P. Buividovich and M. Polikarpov, Entanglement Entropy in Gauge Theories and the 
Holographic Principle for Electric Strings, Phys.Lett. B670 (2008) 141-145, 

[arXiv: 0806.3376], 


26 


[13] Y. Nakagawa, A. Nakamura, S. Motoki, and V. Zakharov, Entanglement Entropy of 
ST/(3) Yang-Mills Theory , PoS LAT2009 (2009) 188, [arXiv: 0911.2596]. 

[14] H. Casini and M. Huerta, A Finite Entanglement Entropy and the C-Theorem, 
Phys.Lett. B600 (2004) 142-150, [hep-th/0405111], 

[15] R. C. Myers and A. Sinha, Seeing a C-Theorem with Holography, Phys. Rev. D82 
(2010) 046006, [arXiv: 1006.1263], 

[16] R. C. Myers and A. Sinha, Holographic C-Theorems in Arbitrary Dimensions, JHEP 
1101 (2011) 125, [arXiv: 1011.5819], 

[17] H. Liu and M. Mezei, A Refinement of Entanglement Entropy and the Number of 
Degrees of Freedom, JHEP 1304 (2013) 162, [arXiv: 1202.2070], 

[18] H. Casini and M. Huerta, On the RG Running of the Entanglement Entropy of a 
Circle, Phys.Rev. D85 (2012) 125016, [arXiv: 1202.5650]. 

[19] H. Casini, M. Huerta, and R. C. Myers, Towards a Derivation of Holographic 
Entanglement Entropy, JHEP 1105 (2011) 036, [arXiv: 1102.0440], 

[20] P. Calabrese, J. Cardy, and E. Tonni, Entanglement Entropy of Two Disjoint Intervals 
in Conformal Field Theory, J.Stat.Mech. 0911 (2009) P11001, [arXiv:0905.2069]. 

[21] P. Calabrese, J. Cardy, and E. Tonni, Entanglement Entropy of Two Disjoint Intervals 
in Conformal Field Theory II, J.Stat.Mech. 1101 (2011) P01021, [arXiv: 1011.5482]. 

[22] M. Headrick, Entanglement Renyi Entropies in Holographic Theories, Phys.Rev. D82 
(2010) 126010, [arXiv: 1006.0047], 

[23] T. Hartman, Entanglement Entropy at Large Central Charge, arXiv : 1303.6955. 

[24] T. Faulkner, The Entanglement Renyi Entropies of Disjoint Intervals in AdS/CFT, 
arXiv:1303.7221. 

[25] V. Rosenhaus and M. Snrolkin, Entanglement Entropy: a Perturbative Calculation, 
arXiv:1403.3733. 

[26] V. Rosenhaus and M. Snrolkin, Entanglement Entropy, Planar Surfaces, and Spectral 
Functions, JHEP 1409 (2014) 119, [arXiv: 1407.2891], 

[27] A. Allais and M. Mezei, Some Results on the Shape Dependence of Entanglement and 
REnyi Entropies, arXiv: 1407.7249. 

[28] A. Lewkowycz and E. Perlmutter, Universality in the Geometric Dependence of Renyi 
Entropy, arXiv: 1407.8171. 

[29] O. Ben-Anri, D. Carnri and J. Sonnenschein, Holographic Entanglement Entropy of 
Multiple Strips, JHEP 1411 (2014) 144, [arXiv: 1409.6305], 

[30] W. Fischler, A. Kundu and S. Kundu, Holographic Mutual Information at Finite 
Temperature, Phys. Rev. D87 (2013) 126012, [arXiv: 1212.4764], 


27 


[31] T. Grover, A. M. Turner, and A. Vishwanath, Entanglement entropy of gapped phases 
and topological order in three dimensions, Physical Review B 84 (2011), no. 19 195120. 

[32] I. R. Klebanov, T. Nishioka, S. S. Pufu, and B. R. Safdi, On Shape Dependence and 
RG Flow of Entanglement Entropy , JHEP 1207 (2012) 001, [arXiv: 1204.4160]. 

[33] B. R. Safdi, Exact and Numerical Results on Entanglement Entropy in 
(5+l)-Dimensional CFT, JHEP 1212 (2012) 005, [arXiv: 1206.5025], 

[34] P. Sabella-Garnier, Mutual Information on the Fuzzy Sphere , arXiv: 1409.7069. 

[35] T. Hirata and T. Takayanagi, AdS/CFT and Strong Subadditivity of Entanglement 
Entropy, JHEP 02 (2007) 042, [hep-th/0608213], 

[36] E. H. Lieb and M. B. Ruskai, Proof of the strong subadditivity of quantum-mechanical 
entropy, Journal of Mathematical Physics 14 (1973), no. 12 1938-1941. 

[37] M. Srednicki, Entropy and Area, Phys.Rev.Lett. 71 (1993) 666-669, [hep-th/9303048], 

[38] N. Shiba, Entanglement Entropy of Two Black Holes and Entanglement Entropic 
Force, Phys.Rev. D83 (2011) 065002, [arXiv: 1011.3760]. 

[39] N. Shiba, Entanglement Entropy of Two Spheres, JHEP 1207 (2012) 100, 

[arXiv: 1201.4865], 

[40] J. Cardy, Some Results on the Mutual Information of Disjoint Regions in Higher 
Dimensions, J.Phys. A46 (2013) 285402, [arXiv: 1304.7985], 

[41] A. Lewkowycz and J. Maldacena, Generalized Gravitational Entropy, JHEP 1308 
(2013) 090, [arXiv: 1304.4926], 

[42] M. Cvetic, G. W. Gibbons, H. Lu, and C. N. Pope, Ricci-Flat Metrics, Harmonic 
Forms and Brane Resolutions, Commun. Math. Phys. 232 (2003) 457-500, 
[hep-th/0012011], 

[43] I. R. Klebanov, S. S. Pufu, and B. R. Safdi, F-Theorem without Supersymmetry, JHEP 
1110 (2011) 038, [arXiv: 1105.4598], 

[44] M. Huerta, Numerical Determination of the Entanglement Entropy for Free Fields in 
the Cylinder, Phys.Lett. B710 (2012) 691-696, [arXiv: 1112.1277], 

[45] H. Casini and M. Huerta, Entanglement and Alpha Entropies for a Massive Scalar 
Field in Two Dimensions, J.Stat.Mech. 0512 (2005) P12012, [cond-mat/0511014]. 

[46] H. Casini and M. Huerta, Entanglement Entropy in Free Quantum Field Theory, J. 
Phys. A42 (2009) 504007, [arXiv:0905.2562], 

[47] N. Drukker and B. Fiol, On the Integrability of Wilson Loops in AdSs x S 5 : Some 
Periodic Ansatze, JHEP 0601 (2006) 056, [hep-th/0506058], 

[48] A. Dekel and T. Klose, Correlation Function of Circular Wilson Loops at Strong 
Coupling, JHEP 1311 (2013) 117, [arXiv: 1309.3203], 


28 


[49] P. Fonda, L. Giomi, A. Salvio, and E. Tonni, On Shape Dependence of Holographic 
Mutual Information in AdSf, arXiv: 1411.3608. 

[50] T. Faulkner, A. Lewkowycz, and J. Maldacena, Quantum Corrections to Holographic 
Entanglement Entropy , JHEP 1311 (2013) 074, [arXiv: 1307.2892]. 

[51] N. Engelhardt and A. C. Wall, Quantum Extremal Surfaces: Holographic Entanglement 
Entropy Beyond the Classical Regime , arXiv: 1408.3203. 

[52] M. Headrick and T. Takayanagi, A Holographic Proof of the Strong Subadditivity of 
Entanglement Entropy , Phys.Rev. D76 (2007) 106013, [arXiv:0704.3719]. 

[53] I. R. Klebanov and S. S. Pufu, M-Branes and Metastable States, JHEP 08 (2011) 035, 
[arXiv: 1006.3587], 

[54] C. A. Agon and H. J. Schnitzer, Holographic Mutual Information at small separations, 
arXiv:1501.0377. 

[55] L. Girardello, M. Petrini, M. Porrati, and A. Zaffaroni, Novel Local CFT and Exact 
Results on Perturbations of N = 4 Super Yang-Mills from AdS Dynamics, JHEP 12 
(1998) 022, [hep-th/9810126] . 

[56] D. Z. Freedman, S. S. Gubser, K. Pilch, and N. P. Warner, Renormalization Group 
Flows from Holography Supersymmetry and a C- Theorem, Adv. Theor. Math. Phys. 3 
(1999) 363-417, [hep-th/9904017]. 

[57] C. P. Herzog and M. Spillane, Tracing Through Scalar Entanglement, Phys.Rev. D87 
(2013) 025012, [arXiv: 1209.6368], 

[58] C. P. Herzog and T. Nishioka, Entanglement Entropy of a Massive Fermion on a 
Torus, JHEP 1303 (2013) 077, [arXiv: 1301.0336], 

[59] J. Cardy and C. P. Herzog, Universal Thermal Corrections to Single Interval 
Entanglement Entropy for Conformal Field Theories, Phys.Rev.Lett. 112 (2014) 

171603, [arXiv: 1403.0578], 

[60] C. P. Herzog, Universal Thermal Corrections to Entanglement Entropy for Conformal 
Field Theories on Spheres, JHEP 1410 (2014) 28, [arXiv: 1407.1358], 

[61] C. P. Herzog and J. Nian, Thermal Corrections to REnyi Entropies for Conformal 
Field Theories, arXiv: 1411.6505. 


29 


